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Abstract 

The generation of electron surface oscillations in overdense plasmas irradiated at normal inci- 
dence by an intense laser pulse is investigated. Two-dimensional (2D) particle-in-cell simulations 
show a transition from a planar, electrostatic oscillation at 2uj, with uj the laser frequency, to a 2D 
electromagnetic oscillation at frequency uj and wavevector k > lo/c. A new electron parametric 
instability, involving the decay of a ID electrostatic oscillation into two surface waves, is introduced 
to explain the basic features of the 2D oscillations. This effect leads to the rippling of the plasma 
surface within a few laser cycles, and is likely to have a strong impact on laser interaction with 
solid targets. 

PACS numbers: 52.38.-r,52.38.Dx,52.65.-y 
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The interaction of sub-picosecond, high-intensity laser pulses with solid targets is of great 
relevance to the generation of bright sources of energetic radiation as well as a test bed for 
Fast Ignitor physics. Since for a solid target the electron density ng ^ ric, where nc = 1.1 x 
10^^ cm~'^/[A/ yum]^ is the cut-off density for laser propagation at the wavelength A, the laser- 
plasma coupling occurs at the target surface over a narrow region with a depth of the order of 
the skin length dp <^ A. The laser force on the plasma has both secular components (leading 
to plasma acceleration, profile steepening and hole boring) and oscillating components; these 
latter drive an oscillatory motion of the "critical" surface where Ug = Uc, that acts as a 
"moving mirror" leading to the appearance of high harmonics in the reflected light [0] . 

Experiments [0] and simulations [0, ^, |^ suggest that either pre-imposed or self-generated 
deformations of the target surface strongly affect laser energy absorption. Evidence for 
small-scale deformations comes from the wide spreading of the reflected radiation observed 
in experiments |p, ^ || at high intensities (~ 10^^ W cm~^) and even for pulse durations as 
short as 35 fs 0. This suggests that surface rippling involves some "fast" mechanism related 
to electron motion rather than Rayleigh- Taylor-like (RT) hydrodynamic instabilities driven 
by the strong target acceleration and occuring on time scales of the ion motion. Seeding 
of density "ripples" by electron instabilities in the underdense plasma region in front of the 
target was evidenced in simulations |^ and investigated theoretically in 0. 

In this Letter, we show with numerical simulations that electron surface oscillations 
(ESOs) may grow for a step density profile (no underdense region present) much faster 
than the typical time scale of ion motion, leading to an oscillatory "rippling" of the critical 
surface. To interpret the simulations we present a model for a new parametric instability, 
based on the "decay" of pumped one-dimensional electrostatic oscillations into two electron 
surface waves. This represents a new nonlinear mechanism of surface wave excitation by 



laser pulses different from previously investigated models fT^- Its potentially strong im- 
pact on the disruption of "moving mirrors" for high harmonics generation (HHG) and the 
production of fast electron jets are discussed. 

We use two-dimensional (2D) particle-in-cell (PIC) simulations to study the dynamics of 
the ESOs with proper spatial and temporal resolution. In particular, in order to evaluate 
the frequency of the ESOs, the complete output of 2D fields was produced eight times for 
each laser cycle. In the simulations reported, the laser pulse is normally incident, impinges 
from the left on the x-axis, and has a wavelength A = 0.25/im [ric = 1.6 x 10^^ cm~^), a 
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uniform spatial profile in the transverse (y) direction and a temporal profile that rises for 
three cycles and then remains constant. The laser pulse is "s-polarized" , i.e. its electric 
field is in the z direction normal to the simulation plane. The plasma has immobile ions, 
initial temperature Tg = 5 keV and a step-like density profile filling the right part of the 
box (x > 0) A numerical box 12A x 4A is taken with a spatial resolution equal to the Debye 
length Xd = \jTe/4:nnoe'^, where Uo is the initial density, and 25 particles per cell are used. 
We will focus on two typical simulations (run 1 and run 2). Introducing the dimensionless 
irradiance Oq = 0.85(/A^/10^^ W/im^/cm^)^/^ in run 1 we take = 1.7 and Uo = Sn^, while 
in run 2 we take Oo = 0.85 and no = Suc- 

The contours of the electron density ne(x, y) in Fig. |l]for run 1 at times t = 8, 10, 12, and 
14 laser cycles from the run start, show the evolution of the surface dynamics over several 
laser cycles. Correspondingly, the space-time contours of ne{x,y = yi,t) at yi = 2. OA and 
?/2 = 1.875A in Fig. ^ show the temporal behavior of the surface oscillations. Initially, the 
surface oscillation is planar, i.e. uniform along y, and has a frequency 2uj, being u = 2ttc/ A 
the laser frequency. It is natural to identify this ID motion as the "moving mirror" driven 
by the longitudinal j x B force at 2uj. In the compression phase, electrons pile up in a 
narrow layer where the peak density is ~ 2no] the electrostatic field E^'^'' (not shown) is 
positive and counteracts the j x B force. In the expansion phase electrons are dragged out 
into vacuum, forming a "cloud" of negative charge with a negative electrostatic field. 

The growth of surface "ripples" can be observed in Fig.[I|. At t = 10, they have small 
wavelengths (~ O.IA), while at t > 12 they evolve into a steady oscillation with wavelength 
Xs ~ 0.5A (Fig.p and frequency ^ uj, which is superimposed to the oscillation at 2uj (Fig. 
This "period doubling" effect is also observed in run 2, for which the density contours at 
t = 16.5, 17, 17.5, and 18.0 are shown in Fig. |^: it is evident that the curvature of the density 
layer is inverted each half laser cycle. In run 2 the oscillation amplitude and the density 
compression is lower than in run 1, and the deformation wavelength is larger (As — 0.75A). 
In the following we will discuss only the long-wavelength structures oscillating at uj and refer 
to them as (2D)ESOs. The 2DESOs are "standing", i.e. not propagating along y. From 
Fig.0 we see that while at yi the amplitude is close to its maximum, at y2 = yi — Xs/4 
there is no evident growth of the oscillation at u, while a weakening of the 2uj oscillation is 
observed. Results (not shown) for p-polarization at normal incidence (i.e, with the electric 
field of the laser pulse in the y direction) are qualitatively similar, the main difference being 
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that the plasma "plumes" extending into vacuum are bent by the laser field in the plane. 

To our knowledge, there are no previous simulations or models describing surface "ripples" 
of the electron density oscillating at the laser frequency for normal incidence. In the "moving 
mirror" motion a superposition of uj and 2uj motions occurs only for oblique incidence and 
p-polarization, because in such a case both the electric and magnetic forces have components 



normal to the surface In ref.|T^ grating- like surface inhomogeneities oscillating at 2uj 
and induced by the magnetic force at oblique incidence were studied. Stationary surface 
ripples in the ion density were observed in simulations |§] at much longer times [|l^ and 
interpreted on the basis of a hydrodynamic model for an underdense, homogeneous, neutral 
plasma which leads (for s-polarization only) to an instability pumped by the secular part of 
the ponderomotive force. Our explanation is that the 2DESOs are generated by a parametric 
decay of the forced ID oscillation with frequency 2uj and transverse wavevector ky = into 
two electron surface waves (ESWs) {u!i,ki) and {u2,k2). Using the matching conditions 
ki + k2 = and Ui + uj2 = 2uj we find that the two overlapping ESWs thus form a standing 
oscillation with frequency uj = ui = 002 and wavevector k = ki = —k2. 

A cold fiuid, non-relativistic 2D model of the parametric excitation of ESWs has 
been developed. Here we report only a brief description of this model while a detailed 
derivation will be published elsewhere |TT]. Firstly, the quiver motion along z is solved 
and a system of 2D Maxwell-Euler equations is obtained where the laser action enters 
via the ponderomotive force. The electron fiuid density and velocity are written as 
He = no{x) + t5nf'^\x,t) + e^5n'~^\x,y,t) and Vg = eV^'^'^\x,t)yi + eV('^)(x, ?/, t), where 
e ~ a1{nc/ne) is a small expansion parameter. The terms at 2uj describe the electrostatic, 
ID "moving mirror" oscillation, that acts as a pump for the instability, and is assumed to 
be unperturbed by the ESWs. The terms at uj are the superposition of two ESWs: 

v(") = e-*"* (v+fcC^'^J^ + w_ke-'^y) /2 + c.c. , (1) 

where \±k = v±fc(x,t) varies slowly in time. To order e^, the coupling between ID and 2D 
modes may be neglected, so that one obtains the usual dispersion relation for "H" surface 
waves propagating along a density discontinuity: 

ec' = u;'{u;l-u;')/{u;l-2u;'), (2) 

where uj^ = ATTUoe"^ /rrif. is the plasma frequency. The evanescence length of the ESWs in the 
plasma is Lsw = (c/'^)!! — 2a;^/ci;p)^/''(l — tu^/cUp)""'^/^. Notice that for electron surface waves 



V ■ E^'^) = —4:TTeSn^^\x,y,t) = 0, and that their maximum frequency is Umax = u'p/V^, so 
that the matching conditions may be satisfied only if < 2uJmax, i-e. Ue > 2nc. 

Inserting the value of the laser frequency in we find that the expected wavelength of 
deformations in run 2 is = 27r//c ~ 0.71A, in good agreement with the simulation. For 
run 1 one finds As ~ 0.87A, quite larger than the numerical result. This is not surprising 
since our expansion procedure is not applicable for the parameters of run 1, where the 
interaction is in the relativistic regime. This may cause, for instance, a lowering of the 



effective plasma frequency by the (time-averaged) relativistic factor ^ ^Jl + al/2, mostly 
due to the relativistic quiver motion along z. By replacing Up by ujp/'jo in (0), for run 1 we 
obtain As ~ 0.55A, much closer to the numerical result. 

By keeping only terms up order in the Euler equation and neglecting feedback effects 
on the ID motion, the Euler equation for the ESW velocity is dtv'^^^ = — eE'^'^-'/me + ea^| 
where a^]_ describes the nonlinear coupling with the ID motion: 

aSvl = -\/i^^)9.v(-) - .i-)a.Fi^-)x + ^5(-)y. (3) 

nieC 

Using this equation and Poynting's theorem the rate of growth of the surface energy for the 
2DESOs was evaluated as F = U^^dtU, where U is the energy density per wavelength of the 
two ESWs and the cycle-averaged variation is 

dtU = Jdx (v('^) ■ (e5n(2'^)E('^) + meUodtv^'^^)) 

-men,V;^2-)9,vl, + n4V;(2-)5,V,y ) + c.c. (4) 

Substituting in the integrand for the expressions of the (unperturbed) ESW fields, finally 
one obtains the growth rate as 

(a - 1)3/2 

F ~ 4:UJal— -T--^ -77^ -7777 (5) 

- 4|[(a- 1)2 + l](a- 2)1/2 W 

where a = rie/ric = uJ^juP. The denominator [a — 4) is actually due to the resonant 
excitation of longitudinal plasmons at uop = 2uj, which makes V^'^^^ very large and invalidates 
our ordering assumptions near resonance. We note that F diverges also for a — > 2; however, 
in this limit the ESW wavelength is very small and thus one expects a strong damping by 
thermal effects neglected in the cold fluid model. 
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Our model is valid at normal incidence for both s- and j9-polarizations (in the sense of 
the 2D simulations), since even if the laser electric field is parallel to the fc- vectors of the 
ESWs there is no resonant coupling of the quiver motion with the ESWs. The transverse 
oscillations are however observed more clearly for s- rather than for p-polarization, since in 
the latter case the motions in the laser and the ESW fields overlap in the polarization plane. 
For oblique incidence and ^^-polarization, in addition to the 2uj motion the laser drives the 
moving mirror at frequency u, which may lead to the pumping of two ESW sidebands around 
the frequency uj/2 [|12|. This latter case is directly relevant to experiments on HHG, where 
the surface rippling observed at high intensities may pose a limit to HHG as an efficient 
XUV source. The normal incidence results already indicate their potential relevance to 
the interpretation of this phenomenon. First, the growth rate is maximum in conditions 
"optimal" for HHG, i.e. when the moving mirror oscillation is driven at high velocities, 
which requires a moderate density "shelf" rather than solid densities. In these conditions, 
at the time of maximum pulse intensity the plasma has a finite gradient at the surface Ue = 
which is in turn steepened by the radiation pressure. Therefore, the plasma profile may look 
rather similar to that assumed in our simulations. In such conditions the ESO instability can 
be much faster than RT instabilities. In fact, even for accelerations of order g ~ 10^'' cm/s^ 
as measured in this regime [|I^, the typical RT growth rate Vrt — V^Rrg — (140 fs)"^ for 
27r/A;/jT ^ A/2 = 0.125/im, is much slower than the rise of the 2DESOs which occurs over a 
few fs in the simulations. 

Finally, the results of run 1 give an indication on the scaling of the growth rate in the 
relativistic regime, not accessible to our analytical model; the growth rate may increase 
strongly for relativistic intensities due to the decrease of the effective plasma frequency, 
producing a stronger rippling of the surface, as consistent with the simulation results and 
experiments that suggest that the instability threshold is close to = 1. 

The ESOs also have a substantial impact on fast electron generation. Fig. ^ shows 
phase space distributions at times t = 8.5 and 13.5 laser cycles for run 1. At early times, the 
momentum distribution is uniform in y, with no accelerated particles in py and most energetic 
electrons having ~ 2meC. At later times, when the 2DESOs have grown, stronger forward 
acceleration occurs near the maxima of the oscillation, showing that most oscillatory energy 
has been transferred to the unstable 2D modes. Correspondingly, strong acceleration in 
Py also occurs. The momentum distribution in px is very regular, with only a minority of 
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electrons "outrunning" the oscillation. This suggests that the generation of fast electrons is 
correlated with the nonlinear evolution and "breaking" of the density oscillations, when the 
amplitude of the latter exceeds the screening length. This may give a spatial "imprint" on 
the transverse structure of the fast electron currents. In previous simulations one may clearly 
observe a spatial correlation between electron jets and "corrugations" at the surface |T^. 
When penetrating into the bulk the jets may either merge or drive current filamentation 



instabilities [|T7|, |T^, thus producing different spatial scales. 

The PIC simulations were performed at the CINECA supercomputing facility (Bologna, 
Italy), sponsored by the INFM supercomputing initiative. 



* Electronic address: macchi@df.unipi.it 
[1] S. V. Bulanov, N. M. Naumova and F. Pegoraro, Phys. Plasmas 1, 745 (1994). 
[2] T. Feurer et al, Phys. Rev. E 56, 4608 (1997). 
[3] S. C. Wilks et al, Phys. Rev. Lett. 69, 1383 (1992). 
[4] H. Ruhl et al, Phys. Rev. Lett. 82, 2095 (1999). 

[5] A. Macchi, H. Ruhl and F. Cornolti, Las. Part. Beams 18, 375 (2000). 

[6] P. A. Norreys et al, Phys. Rev. Lett. 76, 1832 (1996); M. Zepf et al, Phys. Rev. E 58, R5253 
(1998). 

[7] A. Tarasevitch et al, Phys. Rev. E 62, 023816 (2000). 
[8] C. Dietrich et al, 2nd ULIA conference, Pisa (Italy), Sep. 29 - Oct.4, 2000. 
[9] M. G. Cadjan, M. F. Ivanov, and A. V. Ivlev, Phys. Lett. A 222, 324 (1996); Las. Part. Beams 
15, 33 (1997). 

[10] R. Dragila and S. Vukovic, Phys. Rev. Lett. 61, 2759 (1988); E. G. Gamaliy Phys. Rev. E 
48, 516 (1993); S. A. Magnitskii, V. T. Platonenko, and A. V. Tarasishin, AIP Conf. Proc. 
426, 73 (1998), and references therein. 



[11] A. Macchi, F. Cornolti, and F. Pegoraro, preprint physics/0105017 . 

[12] A. Macchi et al, in Superstrong Fields in Plasmas 2, edited by M. Lontano (American Institute 

of Physics, New York, 2002), in press. 
[13] D. von der Linde and K. Rzazewski, Appl. Phys. B 63, 499 (1996). 
[14] L. Plaja, L. Roso, and E. Conejero-Jarque, Astrophys. J. Supp. Ser. 127, 445 (2000). 



7 



FIG. 1: Contours of normalized electron density ne/ric for run 1 (cq = 1.7, ng = Sric) at various 
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target surface is shown. 

FIG. 2: Space-time evolution of ne{x,y = yi,t) at yi = 1.875/A (left) and y2 = 2. OA (right) for 
run 1 (see Figjl|). The position of the Ue = Uc surface is evidenced by a black contour line. 
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FIG. 3: Same as Fig. |l| for run 2 (oo = 0.85, Uq = Stt-c). 



8 



FIG. 4: Phase space projections for run 1 in the {y,Px) (top), iy,Py) (middle), and (x,Px) (bottom) 
planes for t = 8.5 (left) and t = 13.5 (right) laser cycles. 



9 



This figure "Figl_C.gif" is available in "gif" format from: 



http://arXiv.org/ps/physics/0105019v3 



This figure "Fig2_C.gif" is available in "gif" format from: 



http://arXiv.org/ps/physics/0105019v3 



This figure "Fig3_C.gif" is available in "gif" format from: 



http://arXiv.org/ps/physics/0105019v3 



This figure "Fig4_C.gif" is available in "gif" format from: 



http://arXiv.org/ps/physics/0105019v3 



